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ABSTRAGT 

The problem of impulsive heating of dust grains in cold, dense interstellar clouds is revisited the¬ 
oretically, with the aim to better understand leading mechanisms of the explosive desorption of icy 
mantles. It is rigorously shown that if the heating of a reactive medium occurs within a sufficiently 
localized spot (e.g., heating of mantles by cosmic rays), then the subsequent thermal evolution is 
characterized by a single dimensionless number A. This number identifies a bifurcation between two 
distinct regimes: When A exceeds a critical value (threshold), the heat equation exhibits the explosive 
solution, i.e., the thermal (chemical) explosion is triggered. Otherwise, thermal diffusion causes the 
deposited heat to spread over the entire grain - this regime is commonly known as the whole-grain 
heating. The theory allows us to find a critical combination of the physical parameters that govern 
the explosion of icy mantles due to impulsive spot heating. In particular, the calculations suggest that 
heavy cosmic ray species (e.g., iron ions) colliding with dust are able to trigger the explosion. Based 
on the recently calculated local cosmic-ray spectra, the expected rate of the explosive desorption is 
estimated. The efficiency of the desorption, which affects all solid species independent of their binding 
energy, is shown to be comparable with other cosmic-ray desorption mechanisms typically considered 
in the literature. Also, the theory allows us to estimate maximum abundances of reactive species 
that may be stored in the mantles, which provides important constraints on available astrochemical 
models. 

Subject headings: ISM: dust - ISM: clouds - astrochemistry - ISM: cosmic rays 


1. INTRODUCTION 


The earliest stages of star formation occur in cold 
(T ^ 10 K), dense {n{H) > 10“^ cm~^ ), and dark 
(Ay > 10 mag) molecular cloud cores le.g.. iMvers et al.l 
Il987fl . Under such physical conditions, rapid freeze-out 
of molecular species from the gas phase on interstellar 
grains should occur on a timescale of ^ 10®/n(iJ) yrs. 
However, while infrared observations confirm the exis- 
tence of thick ic y mantles on interstellar grains (e.g., 
iGibb et HI 120041 1 , mol ecular species are also observed 
in the dark cold gas dTafalla et al.l 120021: iGaselli et al.l 
I 200 I iTafalla et all 12004 IGaselli et al.l I20m As the 
lifetime of cold mole cular cores is at least ^ 10® yr 
(iBriinken et al.lF2014[ l. a non-thermal desorption mech¬ 
anism is required to maintain the observed gas-phase 
abundances of species. The recent discovery of complex 
organic molecules (lOberg et al.l 120101: IBacmann et al.l 
120121: iGe rnicharo et al.l 120121 1 and deuterated methanol 
(jBizzocchi et al.ll2014ll in the cold gas is further evidence 
for non-t hermal processing and eva poration of cold icy 
mantles (iVasvunin fc Herbstll2013bll . 

Interactions of interstellar grains with cosmic ray 
(CR) particles. X-ray and UV photons, and even 
their mutual collisions cause the grain heating and 
hence stimulate sublimat i on of ice (Id’Hendecourt et al.l 
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heated region may be localized or it may extend over 
the entire grain - these two scenarios are usually 
referred to as “spot heatin g ” and “ whole-grain heating”, 
respectively (iLeger et al.l 119851: iSchutte fc GreenbergI 
119911 : iShen et al.l 12004 iBringa fc Johnsonl 1200411 . Also, 
some exothermic reactions occurring on the grain 
surface (e.g., the formation of molecular hydrogen) may 
result in the local heating and lea d to the chemical 
desorption of weakly bound species (iDulev fc Williamd 
19931: iGarrod et all 120071: iCecchi-Pestellini et al.l 1201^ 
Rawlings et al.ll2013ll . 

One can identify two distinct regimes of desorp¬ 
tion occurring in response to the impulsive grain 
heating: The classical thermal evaporation, and 

the so-called “explosive desorption” triggered by the 
exothermic chemical reaction(s) between free radi¬ 
cals f rozen in th e bulk of ice (Id’Hendecourt et ^ 
19821: ILeger et ahl 119851: ISchutte fc GreenbergI 119911 : 
Shalabiea fc GreenbergI 119941 1. The essential difference 
between the two regimes is that the evaporation of the 
ice mantle (typically limited to the most volatile species) 
is accompanied by the grain cooling, whereas the chemi¬ 
cal reactions (activated by the deposited energy) can lead 
to the runaway temperature growth. As a result, the ex¬ 
plosive desorption may cause the ejection of the entire 
mantle off the grain surface. 

Since the 1980’s, there have been various mechanisms 
proposed to trigger the thermal explosion of icy man¬ 
tles. In particular, these include inelastic collisions 
between the grains, when a certain fraction of their 
kinetic energy is convert e d into heat in the mantle 
(Id’Hendeconrt et'HI 119821: ISchutte fc GreenbergI 119911 : 
iShalabiea fc GreenbergI Il994i) . and the impact of ener- 
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getic particles, su ch as CR and X rays (|Leger et al.lll985l : 
iShen et ^120041) ■ The analysis, however, has been al¬ 
most completely focused on the whole-grain-heating sce¬ 
nario, neglecting the initial thermal spikes emerging in 
a grain (e.g., along the CR paths). To the best of 
our knowledge, the possibility of thermal e xplosion due 
to CR spot heating was only discussed by iLeeer et"^ 
(|1985ll . who concluded that such process is not feasible Q 
In this article we revisit the problem of spot heating 
of interstellar grains. We introduce a concept of the lo¬ 
calized ignition spot and show that the evolution of the 
initial kinetic energy deposited in a reactive medium in 
this case is uniquely described by a single dimensionless 
number A. This concept allows us to calculate a critical 
value of A above which the thermal explosion is triggered 
and, hence, to find a critical combination of the physical 
parameters that govern the explosion of icy mantles due 
to spot heating. We show that the energy deposited by 
iron CR are sufficient to cause such explosions. Further¬ 
more, we demonstrate that the chemical explosion due to 
whole-grain heating is inhibited by efficient sublimation 
cooling. Based on the recent calculations of local CR 
energy spectra, we obtain the minimum expected rate 
of mantle disruption due to impacts of iron CR. Finally, 
the presented theory allows us to estimate the maximum 
abundances of reactive species that may be stored in the 
mantles and, thus, to impose important constraints on 
available astrochemical models. 


2. THEORY 

Consider the situation when a certain amount of ki¬ 
netic energy is “instantaneously” deposited in a reactive 
medium. It is intuitive to expect that the exact form 
of the initial energy distribution must be unimportant 
for its subsequent evolution, provided this distribution 
is sufficiently localized and the energy is rapidly ther- 
malized. Mathematically, the possibility of the thermal 
explosion in this case can be investigated by assuming the 
initial temperature distribution in the form of the delta 
function. Limits of applicability for such approximation 
of the ignition spot are determined from the numerical 
analysis, as discussed below. 

Let us consider the cases where the initial energy is 
concentrated on a plane, along an axis, or in a point. 
For such ignition spots, the problem is characterized by 
the symmetry indices D = 1,2, and 3, respectively, and 
the heat equation describing the temperature distribu- 
tion T(r,t) in a reactive medium has the following form 
(jLandau fc Lifshit3ll987fl : 


at 


f d'^T 
I 5r2 


D-\dT\ 
r dr ) ' 


( 1 ) 


with the initial condition 

T{r,0) = ^6D{r). ( 2 ) 

pc 

Here, 5d{t) is the delta function in D dimensions, qd is 
the initial energy density in the ignition spot, k, p, and c 
are, respectively, the thermal conductivity, mass density. 


^ We note that non-explosiv e desorption due to spot heating has 
been extensively studied ('e.g.. iLeger et anil98fil : IShen et al.ll2004l : 
IBringa Hr. ,lohnsonll2nn4ll . 



Fig. 1. — Evolution of the dimensionless temperature at the cen¬ 
ter of the ignition spot, 9(0, t). Shown are numerical solutions 
of Eqs. 0 and for D = 2, obtained for marginally under- 
critical (dashed line) and over-critical (solid line) values of A. The 
under-critical curve tends asymptotically to the solution for A = 0, 
d(0, r) oc (thin solid line), the over-critical curve approaches 

0(0, r) oc T after the explosive growth. For D = 1 and 3 a similar 
behavior is observed. 



Fig. 2. — Dependence of the numerically calculated explosion 
threshold, Acr, on the normalized size of the ignition spot, w. The 
thresholds, plotted for D = 1,2, and 3, are nearly constant for w < 
1, and rapidly increase at larger w (see Appendix El for details). 
The numerical accuracy for Acr is better than ±5%. 


and specific heat of the medium (treated as incompress¬ 
ible, so c should be taken at constant pressure), is 
the heat of reaction per unit volume and time, and Fla 
is the relevant activation energy in the Arrhenius factor 
(definition of Qr and proper choice of are discussed in 
Sec. 13.11) . We first assume the properties of the medium 
to be independent of the temperature - the cases when 
c or/and k are functions of T are considered later. 

For the analysis of Eqs. (HD and m we normalize the 
temperature by the activation energy, 9 = k^T/E^,. For 
the dimensionless distance ^ = r/r, we choose the scale 
r, which provides unity normalization of Eq. while 
the dimensionless time r = t/t* is determined by the 
timescale t* of thermal diffusion at the distance r*. This 
yields 


r* 



t 


( 3 ) 


so the heat equation is reduced to 

89 8^9 D-189 

— = Ae“^/^ -I-1- 

8t 8^ C dC 


( 4 ) 
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Fig. 3.— Effect of the temperature-dependent specific heat on 
the explosion threshold. The critical numbers Acr, numerically 
calculated for D = 1, 2, and 3, are plotted versus the exponent a. 
determining the temperature dependence c oc T“. The numerical 
accuracy for Acr is better than ±5%. 


and the initial condition - to 


0(e,O)=fe(C). (5) 


Thus, in the dimensionless form the problem is charac¬ 
terized by a single number, 


Qv f qp 

kE^ \pcE^) 


( 6 ) 


The role of A is similar t o tha t of the Frank-Kamenetskii 
number Afk [see Eq. m] which governs the ther¬ 
mal stability of a steady state (|Frank-Kamenetskiilll969l : 
iLandau fc Lifshita Il987fl . The relation between A and 
Afk is discussed in Appendix [XI 
The thermal explosion is triggered when A exceeds a 
certain critical value Acr “ the explosion threshold. From 
the numerical solution of Eqs. o and (O we obtain the 
following thresholds: 


D = l: Acr = 1-45; 
D = 2 : Acr = 9-94; 
D = 3 : Acr = 22.1. 


The bifurcation between the decaying and explosive evo¬ 
lutions is illustrated for D = 2 in Fig. [l] where the 
temperature at the center of the ignition spot, 0(0, r), 
is plotted. For A < Acr the integral effect of thermal 
diffusion is stronger than that of reaction heating, so 
the asymptotic temperature decay is described by the 
fundamental solution of the heat equation in free space 
([Landau fc Lifshit3ll987ll . which yields 0(0, r) oc . 
When A > Acr, thermal diffusion becomes asymptotically 
negligible and the temperature approaches linear growth, 
since the Arrhenius term in Eq. 0 tends to a constant 
(A) for large 0. In Fig. [H the bifurcation occurs at r ^ 3 
(while for D = 1 and 3 it is at r ^ 10 and ~ 1, respec¬ 
tively). We conclude that the explosion develops within 
the physical time of a few t*. 

For the numerical solution, we approximate the initial 
energy distribution by a rectangular function with width 
w (the delta function formally corresponds to the limit 
ic —>■ 0). The obtained dependence Xct{w) is plotted 
in Fig. [21 showing that the explosion thresholds remain 
practically constant for u; ;< 1. Thus, the problem does 
not (practically) depend on the physical size of the ig¬ 
nition spot as long as it is smaller than ~ r*, i.e., the 


initial energy distribntion for such localized spots is well 
represented by the delta function. 

Once the explosion is triggered, the hot reactive zone 
starts expanding away from the ignition spot. As dis¬ 
cussed in Appendix |b 1 the flame front prop agates with 
a constant speed U determined by Eq. (IBID . In Sec. 13.21 
we demonstrate that the magnitude of U is much smaller 
that the typical sound speed in solids. 

The above results can be generalized for the case where 
properties of the medium depend on the temperature. In 
Appendix |C] we show that for a power-law temperature 
dependence of the specific heat, c{T) oc T°‘ (typical for 
solids, see Sec. IXTD . the explosion threshold rapidly de¬ 
creases with the exponent a. Figure |3l illustrating the 
case D = 2, demonstrates that for the linear tempera¬ 
ture dependence the valne of Acr decreases by one order 
of magnitude, and for the quadratic - by two. Note that 
the number A as we ll as t he f ront speed U in this case are 
given by Eqs. (EU and (IC2I) . We also demonstrate that 
the temperature dependence of the thermal diffusivity 
X = k/ pc has relatively weak effect on the results. 

3. IMPLICATION FOR INTERSTELLAR DUST GRAINS 

In this section, the theory presented in Sec.|2]is applied 
to interstellar dust grains, to obtain conditions when im¬ 
pulsive heating by energetic particles is expected to cause 
the thermal explosion of icy mantles. 

The impulsive heating by CR particles, sketched in 
Fig. m has the axial symmetry and is described by the 
solution for D = 2. The initial energy density g2, which 
enters the dimensionless number A in this case, is equal to 
the stopping power of a CR particle. The stopping power 
depends on the particle kinetic energy e (per n ucleon) 
and exhibits a broad maximum at e = (jZieglerl 

II980I1 : for protons, ~ 0.1 MeV and 92 (eh^) 
10“^° J cm“^; for iron ions, ~ 1 MeV/nucleon and 
02(eFe’^’') J cm“^. The heating by X rays is better 

described by the spherically-symmetric solution, D = 3 
(see discussion in Sec. 1321). 

3.1. Properties of icy mantles 

Let us summarize typical physical properties of man¬ 
tles which determine the magnitude of A. 

For many amorphous solids (including ice) the specific 
heat c increases approximately as oc at lower temper¬ 
atures, with 10“^ J cm“^K“^ ^ pc < 0.3 J cm“^K“^ 
for 10 K < T < 50 K; the growth becomes slower at 
higher temperatures , pc ^ 3 J cm~^K ~^ at T ^ 10^ K 
([Zeller fc: PohllllOTlI : iLeeer et'aI1ll985ll . We employ this 
generic dependence for the estimates below. For the 
thermal conducti vity k = pcx we use the diffusivity 
Y ^ cm^s~^ (Id Hen decourt et al.lll982l : iLeger et all 
1 198,51 iSchiitte fc Creenberglll99lll : the latter is approxi- 
mately constant for many amorphous solids at T > 30 K 
([Zeller fc Pohll[l97l!) . No te that y may decre ase with T 
for amorphous water ice (lAndersson fc Su^ 1200211 . but 
this should only have a minor effect on results (see Ap¬ 
pendix 10. 

For the sake of clarity we suppose that among the va¬ 
riety of reactive species (radicals) stored in the mantle, 
there is a pair (A and B) whose exothermic reaction dom¬ 
inates the heat release (the approach can be straightfor¬ 
wardly generalized to multiple reactions). The heat rate 
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Fig. 4.— Schematic representation of the chemical explosion of icy mantle due to cosmic ray impact, (a) A refractory core of a grain 
is covered by a thick icy mantle, where some reactive species (radicals, blue dots) are stored in the bulk. Each collision of a CR particle 
with a grain is accompanied by loss of energy, deposited along the CR path. This creates a hot, narrow cylindrical region whose 
subsequent evolution is governed by the dimensionless number A, given by Eq. |6)l: If A is below a certain critical value, the deposited 
energy is simply redistributed over the grain’s volume (the whole-grain-heating scenario). Otherwise, the thermal explosion is triggered (c) 
- runaway exothermic reactions generate a cylindrical flame front in the mantle, leading to its disruption. 


is then given by (jLee-er et al.lll985li Qr ~ ETiipxtpBNv, 
where E’r ^ 3 eV is the typical energy release per re¬ 
action, </9A.B = Na,b/N are the fractional abundances 
of the species, v ~ 2 x 10^^ s“^ is their characteristic 
vibration frequency, and TV ~ 3 x 10^ ^ cm~^ is the to¬ 
tal number density of mole cules in ice (|Leger et al.lll985l : 
iSchutte fc Greenberelll99lh . Thus, to estimate the mag¬ 
nitude of Qr we need to know the abundance of reactive 
species. Let us elaborate on this point. 

Direct infrared observations of interstellar ices can only 
supply us with the abundances of major ice constituents 
which are in general not reactive under cold ISM condi¬ 
tions, with the exception of CO ice. As such, we have to 
rely on astrochemical modeling when estimating abun¬ 
dances of reactive species in a typical interstellar ice. 
Early astrochemical models did not have a distinction 
between rea ctive surface and more inert bulk of a thick 
icy m antle (|Hasegawa et al.l 119921 : iHasegawa &: Herbs^ 
1199311 . Therefore, in these models all species adsorbed 
on a grain surface participate in efficient “surface” chem¬ 
istry and the resulting fraction of radicals stored in the 
mantle is very low. However, in a number of more re¬ 
cent studies, several important effects were recognized 
that favor larger amounts of radicals to be stored in 
the interstellar ice. First, icy mantles in dark clouds 
are likely to be thick and consisting of several hundreds 
of mo nolayers (see, e.g.. Sec. 4.2 of iCaselli fc Ceccarellil 
l2012tl . For this reason, reactive species in the inner lay¬ 
ers of the ice may be quickly covered by new accreting 
species during the ice formation, and become excluded 
from the rapid surface chemistry. Reactive species be¬ 
come froz en into water ice an d thus survive and accumu¬ 
late le.g.. lTaQuet et al1l2012tl . Moreover, it is likely that 
icy mantles are exposed to the UV photons even in dark 
clouds. Photons can penetrate the entire mantle and dis¬ 
sociate stable molecule s in the ice (|Cruz-Diaz et al.ll20l4 
Chang fc Herbstir2014ll . thus producing radicals (|Garro(l 


2013a|i . Finally, the amount of radicals in the ice may be 


affected by the internal ice structure: theoretical studies 
show that the p orous structure of ic e favors accumulation 
of the radicals (|Taauet et al.ll201^ . However, some au¬ 
thors sh ow that interst ellar ices are rather compact than 
porous (iGarro dl l2ni3bD . 

To obtain quantitative estimates of the fraction of re¬ 
active species stored in the ice in a dark cloud, we sim¬ 
ulate the formation of the icy mantle during the con¬ 


traction of a diffuse cloud into a dense core using our 
MONAG O code and a simp l e evolu tionary model pre¬ 
sented in iVasvunin fc HerbstI (|2ni3al l. Briefly, in the evo¬ 
lutionary model, the temperature linearly decreases with 
time from 20 K to 10 K, and the gas density increases 
from 10^ cm“^ in the beginning to 10® cm“® at the end 
of the contraction. Visual extinction Ay increases self- 
consistently with the density, from Ay = 3 to Ay > 10. 
The MONAGO cod e has be en updated in comparison to 
iVasvunin fc HerbstI (l2013al) . and now it includes chem¬ 
istry in the bulk, due to ice photoprocessing and intra¬ 
mantle diffusion of species (details will be described in a 
future paper). 

Mobility of species in the bulk of ice is likely to 
be significantly lower than on the surface, due to the 
larger number of neig hboring species that bonded to each 
other (lGarrodll2013ari . In our microscopic formalism, this 
means higher diffusion energy f or a spe c ies in the bulk 
than on the surface. Following iGarrodI (j2013a[ l. we set 
the diffusion energy of species in the bulk to be two times 
the respective surface diffusion energy. The latter, in 
turn, is usually taken as a fraction of the sublimation 
enth alpy (typically, their ratio varies from 0.3 to 0 .8, see, 
e.g.. IHasegawa et all 119921 : iRuffle fc HerbstI 1200011 . here 
we chose the value of 0.5 in agreeme nt with the best-fit 
model bv IVasvunin fc HerbstI (j2013ari . As such, the diffu¬ 
sion energy is ~ 1150 K for GO molecules, so GO as well 
as other abundant reactive species (with higher diffusion 
energies) that are produced during ice photoprocessing 
and entrapment of accreting material can effectively ac¬ 
cumulate in the bulk of ice. Thus, we shall consider the 
GO diffusion energy as the relevant activation energy for 
the Arrhenius factor in Eq. o, i.e., Aa/fcB = 1150 K. 

In Fig. [5l the fractional abundances (p of the most 
abundant reactive species in the ice are plotted ver¬ 
sus time. GO is mainly accreted from the gas phase. 
Some of it undergoes hydrogenation and ultimately con¬ 
verts to methanol and other saturated species, but sig¬ 
nificant fraction of GO molecules get buried in the icy 
mantle in a pristine form. The next most abundant 
species is OH, which is mainly produced via dissocia- 
tion of water by photons and GR protons (according to 
[Andersson fc van Dishoe^ (l2008f) . only a fraction of the 
dissociation products recombine back to H 2 O). Note that 
the abundance of OH in our model is about two orders 
of magnitude lower than in other models of multilayer 
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Fig. 5. — Evolution of the fractional abundances of the major 
reactive species (with respect to the total number of molecules 
in the ice). Results are from the numerical modeling with the 
modified MONACO code. 

ice fe.g.. iTaauet et alll2012D . Presumably, this is due to 
the fact that we take into account efhcient recombination 
of OH with free H atoms that are generated in the bulk 
of ice and perform a random walk before reaching the ice 
surface. Finally, a certain fraction of HCO is produced 
in the bulk, mainly in dissociation of methanol by cosmic 
ray protons. We see that the abundances of CO and OH 
reach the values of ^pco ^ 10“^ and (^oh 3 x 10“^ at 
later stages of the contraction. We employ these charac¬ 
teristic values for the estimates below. 

3.2. Explosion due to spot heating 

In Sec. [2] we pointed out that the presented theory can 
be used as long as the physical size of the ignition spot 
does not exceed ^ r*. By substituting typical parame¬ 
ters (listed above, with pc = 0.3 J cm“^K“^) in Eq. ([3]) 
for D = 2, we obtain r* ^ 3 x 10“® cm for the heating by 
iron CR. This value is substantially larger than the di¬ 
ameter of the cylindrical volume where CR deposit their 
energy (< 100 A, iLeeer et'all 1198511 and, at the same 
time, is smaller than the size of large grains dominating 
the interstellar dust mass 10“^ cm). Furthermore, 
the cylindrical explosion develops during the time of the 
order of 3t* ^ 3 x 10“® s, which is much longer than the 
time during whi ch the deposited C R energy is thermal- 
ized (< 10“^^ s. iLeger et al.lIlQS'SIl . Thus, the theory is 
indeed applicable to study the reaction of large grains on 
impulsive heating by heavy CR species (assuming dust 
properties that are typically used in astrochemical mod¬ 
eling, see Sec. [XT]) . 

Let us estimate the magnitude of A for individual col¬ 
lisions with iron CR. First, we assume a constant c 
(and k) for icy mantles. By substituting in Eq. ([6]) 
q 2 ^ 10“® J cm“^ and pc = 0.3 J cm“^K“^, and set¬ 
ting Fr = 3 eV and (ficofOH = 3 x 10“^ for the reaction 
between CO and OH, we obtain A ^ 30, which exceeds 
Acr ~ 10 for D = 2. Hence, iron ions with the energy 
corresponding to the maximum of the stopping power are 
able to trigger the explosion^ whereas CR protons with 
q 2 ^ 10“^° J cm“^ remain under-critical, since A oc 92 - 

^ The explosion threshold could also be evaluated by directly 
comparing the rates of chemical reaction and thermal diffusion in 
Eq. JTJ. However, as shown in Appendix]^ this comparison should 


Remarkably, when the temperature dependence of the 
specific heat is taken into account, the resulting ratio 
A/Acr becomes even larger, i.e., the explosion condition 
is relaxed in comparison with the constant-c case: For 
c{T) oc T“, we employ the results of Appendix (Cl and cal¬ 
culate the enthalpy scale He with pc{E) ~ 3 J cm“^K“^ 
and 1 < a < 2; using the dependence Acr(ct) plotted in 
Fig. [21 and substituting He in Eq. (ED we obtain A/Acr 
in the range between ^ 3 and ^ 30 for iron CR. Thus, 
even if some numbers used above for estimating A would 
be somewhat less favorable (e.g., if ipA^^B 3 x 10“®), 
iron CR should still lead to the explosion. 

The flame front, generated in the mantle by the explo¬ 
sion, propagates with the speed U given by Eqs. m or 
(IC2I) . From this we obtain U ^ 10'^ cm s“^, which is 
more than an order of mag nitude lower than the typical 
sound speed in ice (see, e.g.. lVogt et al.ll2008ll . The cross¬ 
ing time in a grain of the radius a is ^ a/U ^ 10“® s for 
a ^ 10“® cm, so one would expect a practically instant 
evaporation of the whole mantle. Yet one should keep in 
mind that the flame front exerts enormous stress - the 
thermal pressure ^ NT substantially exceeds GPa-level, 
while the tensi le strength of ice is less than one MPa (e.g., 
IPetrovidl2003il . This may lead to mechanical disruption 
of the mantle before it is completely evaporated. 

One can estimate the rate of mantle disruption due 
to the thermal explosions, 1/tdis, which is determined 
by the local energy spectrum of iron CR. We assume 
a constant abundance of i ron ions of (fipc 10“^ (rel¬ 
ative to protons, see, e.g., ILeger et'D] 119851 : IShen et p] 
200411 and empl o y the local proton spectrum Jh (e) from 
Padovani et al.l ll2009f l . where e is the energy per nu¬ 
cleon. The disruption rate is equal to the product of 
the grain cross section and the CR flux contributing to 
the explosion. The minimum value of the latter can 
be roughly estimated at0 47r())FeeF(f'’^>/H(£F«f’')i where 
~ 1 MeV/nucleon correspo nds to the m aximum 
of stopping power for iron ions ([Ziegler] Il980f) . Even 
for dense clouds (with the column density of molecu¬ 
lar hydrogen of ^ 3 x 10^^ cm~®, where th e spectrum 
is strongly attenuated. iPadovani et all 1200911 . we obtain 
that the disruption rate for large grains (a ^ 10“® cm) 
is not lower than 

1/tdis ~ {2na)^<PF,e^rM£F:i - 10“® 

Furthermore, supposing the entire mantle evaporated 
upon disruption, we can also estimate the minimum 
desorption rate of molecules into the gas phase. For a 
mantle with thickness Aa, the desorption rate of species 
A is ~ Aira'^AaNpA/tdis oc a® (assuming Aa oc a). We 
see that the explosive desorption is heavily dominated 
by large grains, with the desorption rate of the order 
of 3 X 10“^ molecules grain“^s“^ for CO molecules. 
This value is comparable to the desorption rates due to 

be performed at the “optimum moment” (short before the bifur¬ 
cation in Fi g. [Til. Presumably, this latter point was not taken into 
account by IL^er et al.l Ill985lj , who concluded that the thermal 
explosion due to spot heating is unlikely. 

^ To obtain the CR flux contributing to the explosion, one should 
integrate JH(e() over the range of energies (around where 

q 2 {£) exceeds the critical value, i.e., where A(g 2 ) 10- Given the 

uncer tainties in the local spectrum in this range JPadovani et al.l 
[200^ . only the lower bound of the flux can be reasonably esti¬ 
mated. 
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a combination of other mechanisms (|Shen et al.l 12004 
iHerbst fc CuDPenll2006f) . such as explosion due to whole- 
grain heating (see next section for its critical discussion), 
and evaporation due to whole-grain and spot heating. 
We note that t he calculations of the desorption rate re¬ 
porte d earlier (iLeger et m il985t iHartauist fc William^ 
19901: iHaseeawa fc HerbstI Il993t iShen et alJ 12004 
Bringa fc Johnsonll2004 IHerbst fc CuDDenll2006li do not 
take into account attenuation of the local CR spectrum, 
which is included in our analysis. 

It is noteworthy that we can practically exclude other 
sorts of energetic particles (e.g., X rays or UV photons) 
as possible causes of explosion due to spot heating. To 
demonstrate this, let us consider X rays as the most ener¬ 
getic species among such particles: The maximum energy 
which can be deposited in a grain by an X-ray photon is 
limited by the condition that the stopping range of elec¬ 
trons produced by the photon is smaller than the grain 
size; for a ^ 10“° cm we get t he upper energy limit of 
the order of a few keV (see, e.g. . ILeger et al.l[T985ll . Since 
energetic electrons lose most of the energy at the end of 
their paths, the spherically-symmetric solution is more 
appropriate to describe the problem in this case. For 
D — 3 (and otherwise the same parameters as above), 
from Eq. ([S]) we obtain that the minimum ignition en¬ 
ergy to satisfy the condition A > Acr — 22 is ^ 10° eV, 
which exceeds the maximum deposited energy by about 
two orders of magnitude. 

Finally, the presented theory allows us to impose im¬ 
portant constraints on the fractional abundance of reac¬ 
tive species in icy mantles, and thus to discriminate be¬ 
tween different astrochemical models. In particular, one 
can estimate the upper limit of the abundance of radi¬ 
cals which can be stored in a mantle: For example, some 
models predict that at later stages of the cloud evolution, 
the product of the relative abundances of such radicals 
may be as high as ^ 3 x 10~° (or even higher, see, 

e.g., Schutte fc Greenberdll991l: iShalabiea fc Greenberg! 

11994 iTaouet et al.ll2012l : iChang fc Herbstll20l4 l. Since 
A oc its value for iron CR would then be about 

two orders of magnitude larger than Acr, so the obtained 
abundances could already be marginally sufficient to sat¬ 
isfy the explosion condition for CR protons. However, 
the latter are ^ 10^ more abundant than iron CR, so the 
very possibility of mantle explosion due to impacts of 
CR protons would imply unrealistically high disruption 
rates, of ^ 10“° yr“° or even larger (these exceed the 
freeze-out rates at typical molecular cloud densities, i.e., 
icy mantles simply would not have time to grow). Hence, 
such high abundances of radicals can be ruled out based 
on the explosion theory. 

3.3. On the whole-grain heating 

If the stopping power of an energetic particle collid¬ 
ing with a grain is too low (under-critical), the de¬ 
posited energy is rapidly redistributed over the whole 
grain. Even though the overall temperature increase 
in this case could be only a few tens of degrees, this 
leads to an exponential amplification of chemical heating 
in the entire reactive volume of the grain, with impor¬ 
tant consequences for surface chemistry and the chem¬ 
ical composition of icy mantles. The thermal stability 
in this regime is determined by the global balance be¬ 
tween the volume heating and the surface cooling due 


to thermal radiation and sublimation (fLe s er et al.iri985l : 
iSchutte fc Greenberg|ll99 ll iCuDoen et ahl 200611 . There- 
fore, it has been usually argued that there is a criti¬ 
cal temperature of the whole-grain h eating, above which 
the e x plosion mus t be tr iggered (Id’Hendecourt et al l 
1982t ILeger et ahl 119851: ISchutte fc Greenberd 119911 : 
Shalabiea fc Greenberg! 11994 IShen et al.ll2004l . As we 
pointed out in the introduction, other mechanisms of the 
whole-grain heating, e.g., due to inelastic grain-grain col¬ 
lisions, have also been suggested as a possible cause of 
the explosion. 

Let us consider the global thermal balance for a re¬ 
active spherical grain. The steady-state temperature 
distribution inside the grain is almost homogeneous, 
so the heating power Pheat is approximately the prod¬ 
uct of and the reactive (mantle) volume 

47 ra^Aa. The surface cooling at te mperatures above 

25 K is dominated by su blimation (ILeger et al.l 119851 : 
ISchutte fc Greenberg|ll99lll . The resulting cooling power 
Pcooi is the product of the area 47ra^ and the cooling rate 
AiLsub(27rmfcBP)“°/^Poe“^^®“‘’^^®^, where AH^nh and 
m are the sublimation enthalpy and the mass of evap¬ 
orating molecules, respectively, and pn is the pre-factor 
for th e saturated vapor pressure (|Legeiill9^lLdger et al.l 
1198,^ . By substituting the heat rate Qr = E^lpa<^bNv 
for reactive species A and B, we obtain the heating-to- 
cooling power ratio, 

PrA^AQ:/V27rmfcBr 

Pcooi ^ ^ APsubPO 

which must exceed unity for the temperature to increase 
with time. For estimates we set P = 30 K which ensures 
that, irrespective of poorly known emission efficiency of 
grains, the radiative cooling is negligible (ILeger et ahl 
119851 : ISchutte fc Greenberd 119911 : IShen et al.l 120041 . By 
adopting Pa/^B — APgub/fcB — 1150 K for CO 
molecules, Aa = 2 x 10“° cm for the mantle thickness, 
and po — 10°^ dyne cm“^ for the saturated CO-vapor 
pressure (|Leger et al.lll985l l. we get Pheat/Pcooi 10“^ 
(for the reaction between CO and OH). 

We see that the global grain cooling is much more ef¬ 
ficient than the heating. The temperature could only 
increase with time if the sublimation enthalpy APgub 
would be substantially larger than the activation en¬ 
ergy Pa, say by several hundreds of K. However, as was 
pointed out in Sec. I3T1 these two values are estimated 
to be about the sameQ and therefore it is rather unlikely 
that the whole-grain heating could trigger the thermal 
explosion. 

Under-critical energetic particles colliding with a grain 
may nevertheless stimula te reactions between radicals 
stored in the mantle ('e.g.. lReboussin et al.ll20l4l . To es¬ 
timate this effect (assuming the whole-grain heating), we 
compare the characteristic time scales of the sublimation 
cooling and the chemical reactions. The time to burn the 
characteristic fraction Cp = ^paPb of the major reactive 
species at a given temperature is phem ((p^)“°e'®'*/^'^^. 

Recent studies JGhesauiere et ^I201f]li suggest that the diffu¬ 
sion energy of CO molecules in the bulk ice (which determines the 
magnitude of Ea) might significantly exceed the value of lAHg^ah- 
If so, the ratio Pieat/Pool would be even smaller than estimated 
above. 































































7 


while the cooling time is tcooi E^ep/Pcooi, where E'dep 
is the total energy deposited in a grain. By substituting 
Edep — 2q2a < 10"* eV for CR protons and using other¬ 
wise the same parameters as above (and also taking into 
account the crossover to the radiative cooling at lower 
temperatures), we obtain tcooi/ichem < 3 x 10“®. 

We conclude that the chemical reactions stimulated 
by the whole-grain heating, due to collisions with under- 
critical particles, are many orders of magnitude slower 
than the cooling. Even though the collisions with (under- 
critical) CR protons are ^ ^ 10^ more frequent 

than with (over-critical) iron ions, such reactions are 
not expected to noticeably affect the chemical compo¬ 
sition of the ice mantle at a timescale of the explosive 
disruption (although the abundances of trace species, 
such as complex organ ic molecules, may be changed, e.g., 
IReboussin et al.l [20141) . It must be stressed, however, 
that the above estimates completely neglect effects of 
the local thermal spikes generated in the mantle by CR 
protons, at a timescale of thermal diffusion. The ques¬ 
tion of whether the integral effect of the heterogeneous 
chemistry stimulated by such heating is more profound 
than that due to whole-grain heating requires a separate 
careful study. 


4. CONCLUSIONS 

The main result of this article is that we have iden¬ 
tified the regime of localized spot heating of a reactive 
medium, and developed a rigorous theory describing the 
thermal evolution in this case. The problem is character¬ 
ized by a single dimensionless number A which depends 
on the deposited energy and properties of the medium. 
The theory allows us to determine the universal explo¬ 
sion threshold and accurately describe impulsive heating 
of icy mantles by energetic particles. 

A collision with an over-critical energetic particle 
(when A exceeds the threshold) leads to the thermal ex¬ 
plosion which, in turn, generates the flame front propa¬ 
gating in the mantle and leading to its disruption. We 
showed that heavy CR species, such as iron ions, are 
able to trigger the explosion, while the stopping power 
of the most abundant CR protons is insufficient for that 
(since the stopping power is roughly proportional to the 
squared atomic number of CR ions). Also, we practically 


ruled out other energetic species, e.g., X rays, as possible 
causes of explosion due to impulsive heating. 

It is important to stress that the question of how ex¬ 
actly the disruption occurs - whether the mantle is com¬ 
pletely evaporated due to thermal explosion, or a part 
of it is ejected off the grain in a form of tiny ice pieces 
- remains unclear. Thus, the possibility of partial me¬ 
chanical disruption of the mantle leads to a conclusion 
that interstellar medium may contain solid nanoparticles 
of predominantly water ice. 

Interestingly, the existence of a well-defined explosion 
threshold allows us to estimate also the upper limit of the 
abundance of radicals that may be stored in the mantle: 
For the assumed dust and CR properties, the product of 
the fractional abundances of two major radicals cannot 
exceed the value of ^ 3 x 10“^, to avoid unrealistically 
large desorption rates. Thus, the presented theory en¬ 
ables us to put constraints on astrochemical models. 

When A is below the threshold, the deposited heat is 
quickly redistributed over the entire grain volume, i.e., 
the whole-grain heating scenario is realized. The rates of 
reactions between radicals frozen in the mantle exponen¬ 
tially depend on the temperature, so even a slight tem¬ 
perature increase can dramatically accelerate the release 
of chemical energy in the reactive volume - for this rea¬ 
son, the whole-grain heating has been considered so far 
as the prime possible cause of thermal explosion. How¬ 
ever, we have demonstrated that the explosion is unlikely 
in this case, since the cooling from the grain surface (due 
to sublimation of volatile species) turns out to be very 
efficient. 

The non-explosive chemical processes, induced in the 
mantle by under-critical impulsive heating, represent an¬ 
other very important phenomenon which needs to be fur¬ 
ther investigated. We considered the whole-grain heat¬ 
ing due to CR protons, and demonstrated that in this 
case the abundance of the major reactive species is not 
expected to noticeably change at a timescale of the ex¬ 
plosive disruption (caused by heavy CR species). How¬ 
ever, chemical reactions depend on the local temperature 
and therefore evolve much faster during short transient 
heating events, within small volumes where under-critical 
particles deposit their energy. Careful analysis of such 
heterogeneous chemistry (as opposed to the chemistry 
due to whole-grain heating), and the evaluation of its 
integral effect will be reported in a future paper. 


APPENDIX 
APPENDIX A 

RELATION BETWEEN UNSTEADY AND STEADY PROBLEMS 


Consider a reactive medium which has a characteristic size vq and temperature Tq at the boundary, and assume that 
there is a thermal equilibrium. The stability of su ch steady state is determined by the Frank-Kamenetskii number 
(|Frank-Kamenets^li969l iLandau fc Lifshit3ll987t) . 


Afk 


nksEQ 


(Al) 


which is the ratio of the time scales of thermal diffusion to chemical reaction: When Afk exceeds a certain threshold, 
the diffusive loss cannot compensate for temperature increase due to ongoing reaction and the steady state becomes 
unstable, i.e. , the thermal explosion i s triggered. The thresholds for U = 1,2, and 3 are AFK,cr = 0.88,2, and 3.32, 
respectively (|Frank-Kamenetskiilll969l ). 

For the unsteady problem studied in this paper, the number A plays a role of Afk- To understand their relation, 
let us calculate the “momentary” value of Afk for the unsteady process: The relevant scale for Tq would be the 
temperature T(0,t) at the center of the ignition spot, while for Tq one should substitute the squared diffusion length 
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{qo/P cTqY/^. Then, by employing Eq. ([S]) we get 


Afk/A - 

where 0o(O = ^b 7’(0, t)/Ea. We see that Afk/A is the sole function of 0o and, thus, of t. It attains maximum at 
00 = 2(1 + 1/!?), where Afk/A ^ 1, which identifies the “optimum moment” to trigger the explosion (provided A > Acr)- 
Physically, the optimum comes out because the size of the reactive zone is too small in the beginning (i.e., the time 
scale of thermal diffusion is short), while at later times the reaction becomes exponentially slow. 

The derived relation allows us to obtain the dependence of Acr on the size of the ignition spot w. Using the relation 
Oqw^ ^ 1 (where w is in units of r*), we get the scaling Xcr{w) ^ which provides excellent fit to 

the curves in Fig. [2] at w > 2. 

Thus, unlike the case of localized ignition spot, the unsteady problem for ru > 1 is no longer characterized by a single 
dimensio nless number. A similar p robl em of thermal expl osion of large “hot spots” has been studied numerically in the 
1960’s by iMerzhanov et al.l (Il963fl and iMerzhanovI (|1966ll who showed that, for a given initial size fq and temperature 
To of the spot, the explosion threshold Afk.cf has a logarithmic dependence on Tq. 


APPENDIX B 
FLAME FRONT 

The explosion generates the flame front propagating away from the ignition spot. At sufficiently large times, when the 
front coordinate ^ is much larger than the front thickness, the last term on the rhs of Eq. @ becomes asymptotically 
negligible (i.e., the front curvature is no longer important). Then, by separating the reactive (9 > 9tr) and inert 
(0 < 0tr) zones of the front (|Frank-Kamenetskiil 119691 : iLandau fc Lifsliit^ 1198711 . we can approximately describe the 
temperature profile by the following equation: 


0 > 0tr : 


0 < 0tr : 


36 _ ^20 


where 0tr ~ 1 is the fitting parameter [to be determined from numerical solution of Eqs. (|1]) and ([S])]. We search 
the solution in the form 0 (^,t) = 0(s) with s = ^ — ut, which yields 0(s) = Aie““® — {X/u)s + A 2 for 0 > 0tr and 
0(s) = ^ 36 ““® for 0 < 0tr. By setting 0(0) = 0tr and taking into account that in the reactive zone 0(s) cannot grow 
faster than linearly, we obtain Ai = 0 ; constants A 2 and A 3 are determined from continuity of 0 and 39/3^ at s = 0 . 
We get u = y^A/0tr and 


S<0: 0(s) = -\/A0tr S + 0tr, 

s> 0 : 0(s) = 0trexp , 

the numerical fit yields 0tr — 1-3. In physical units, the front speed. 


U = 




, 0tr(pc)2Ea’ 

is determined by the reactive and transport properties of the medium. 


(Bl) 


APPENDIX C 

THERMAL EXPLOSION WHEN C OR x ARE FUNCTIONS OF T 

Let us consider the case when the specific heat is a functio n of temperature, c(T), w hile the thermal diffusivity 
X = k/ pc is first assumed to be constant. It is convenient ([Landau fc Lifshit3 1198711 to introduce the enthalpy 
H = p J cdT = F{T), noting that F{T) is a single-valued (monotonously increasing) function. Then Eqs. (P) and (p) 
can be rewritten in the following identical form for F [: 


3H 





D-13H\ 
r 3r ) ' 


i!(r,0) = qoSoir), 


where T = F ^{F[) is the inverse function. We introduce the enthalpy scale, FIe = T(Ta), and conclude that the 
problem can be reduced to the dimensionless form of Eqs. (m and ([ 5 ]), where (apart from the Arrhenius term) H/He 
should be substituted for 0; in the Arrhenius term, 0 should be replaced with F~^{H/He), and 


A = 


Qr 

xHe 


2/0 


Qd 

He 


(Cl) 
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Correspondingly, the speed of the flame front is given by 


U = 



(C2) 


For example, for c oc with the exp onent a > 0, we get HIHe = whe re He = pc(£’a)F^a/(l + o;); for a 

constant specific heat (a = 0) Eqs. ()C1[1 and (|C2p are reduced to Eqs. ^ and (jBip . respectively. Figure [3] shows that 
the explosion threshold decreases dramatically with a. 

Also, we analyzed the effect of the temperature-dependent thermal diffusivity x- this case the diffusion term in 
the heat equation becomes nonlinear. From the numerical solution with x ck we obtained the dependencies Acr(—/5) 
that are qualitatively similar to Acr(Qi) shown in Fig.[3](i.e., Acr monotonically increases with /3). However, the relative 
variation of Acr with /3 turns out to be several times smaller than with a, i.e., the effect of x(^) on the explosion 
threshold is substantially weaker than that of c(T). 
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